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Abstract 

We present a new truncation scheme for the Schwinger-Dyson equations of QCD that respects 
gauge invariance at any level of the dressed loop expansion. When apphed to the gluon self- 
energy, it allows for its non-perturbative treatment without compromising the transversality of the 
solution, even when entire sets of diagrams (most notably the ghost loops) are omitted, or treated 
perturbatively. 
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Introduction - The quantitative understanding of the non-perturbative properties of 
Quantum Chromo dynamics (QCD) |l| constitutes still one of the most challenging prob- 
lems in particle physics. The basic building blocks of this theory are the Green's (correla- 
tion) functions of the fundamental degrees of freedom, gluons, quarks, and ghosts. Their 
non-perturbative structure is at the center stage of extensive research that could furnish 
invaluable clues for deciphering the infrared dynamics of QCD. 

Lattice simulations are indispensable in this quest, since they capture all the non- 
perturbative information of the theory. It has become clear by now that the lattice simu- 
lations yield an infrared finite gluon propagator in the Landau gauge. This rather charac- 
teristic behavior has recently been firmly established using lattices with large volumes; in 
addition, the non-perturbative ghost propagator in the same gauge diverges, at a rate that 
deviates only mildly from the tree- level expectation These clean lattice results consti- 
tute a serious challenge for the QCD theorists: obtaining the same results from the theory 
formulated in the continuum is bound to expose a fundamental dynamical mechanism at 
work. 

In the continuous formulation the dynamics of all Green's functions are determined by 
an infinite system of coupled non-linear integral equations known as Schwinger- Dyson equa- 
tions (SDE) [3]. These equations are inherently non-perturbative and can be used to address 
problems related to e.g., chiral symmetry breaking, dynamical mass generation, and forma- 
tion of bound states. Since this system involves an infinite hierarchy of equations, in practice 
one is severely limited in their use, and the need for a self-consistent truncation scheme is 

n 

evident Devising such a scheme, however, is very challenging, especially in the context 
of non-abelian gauge theories, like QCD [5]. The central problem stems from the fact that 
the SDEs are built out of unphysical off-shell Green's functions; thus, the extraction of re- 
liable physical information depends crucially on delicate all-order cancellations, which may 
be inadvertently distorted in the process of the truncation. 

The situation may best exemplified with the SDE of the gluon propagator AQ,^(g). In 
the Feynman gauge. 



A„^(g) = -i 



gap - —j^ A(g ) + —r 



(1) 



H J Q 

where Iiap{<i) = {dap — <la<ii3 / <f')^{,'f') is the gluon self-energy and A~^(g^) = + iYi{q^). 
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FIG. 1: The PT procedure to construct the new SDE of the gluon propagator. External legs ending 
in a gray circle represents background gluons.The corresponding Feynman rules can be found in 

Q. 

The conventional SDE for U^^ reads 
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na/3(g) = ^{ai)af3, (2) 
i=l 

where the diagrams (oj) are shown in Fig{T^. Since the self-energy enters in the latter 
diagrams (white blobs in the same figure), Eq.([2]) constitutes a dynamical equation that can 
in principle determine Ila/j. Due to general arguments based on the Becchi-Rouet-Stora- 
Tyutin (BRST) symmetry Uaisiq) is transverse, i.e. g°TIo/3(g) = 0. Notice, however, 
that enforcing this fundamental property on the rhs of Eq. ([2]) , i.e., through the contraction of 
individual graphs by g", is far from trivial, essentially due to the complicated Slavnov- Taylor 
identities (STI) satisfied by the fully-dressed vertices. As a result, the SDE of Fig{T^ cannot 
be truncated without compromising the transversality of Uaisiq). For example, keeping only 
graphs (ai) and (02) is not correct even at one loop. Adding (03) is still not sufficient for a 
SDE analysis, because (beyond one-loop) ^"[(cii) + (02) -l- (03)]^^ 7^ 0. 

In this letter we present a new truncation scheme for the SDE of (quarkless) QCD that 
respects gauge invariance at any level of the dressed loop expansion. This becomes possible 
due to the drastic modifications implemented to the building blocks of the SD series, i.e. 
the off-shell Green's functions, following the field-theoretic method known as pinch tech- 
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nique (PT) The PT is a well-defined algorithm that exploits systematically the BRST 
symmetry in order to construct new Green's functions endowed with very special properties. 
Most importantly, they satisfy abelian, Ward identities (WI) instead of the usual STIs, have 
correct analytic properties and displays only physical thresholds js]. 

The PT rearrangement gives rise dynamically to a new SD series analogous to the one 
in Eq.([2]), with the following characteristics: the graphs appearing on the rhs are made out 
of new vertices (Fig{T]i), but contain the conventional self-energy IIq/j as before. These new 
vertices correspond precisely to the Feynman rules of the Background Field Method (BFM) 
in the Feynman gauge, i.e., it is as if the external gluon had been converted dynamically into 
a background gluon. The Ihs, in addition to the term Uapiq) already there, contains addi- 
tional terms, also proportional to Uapiq), which are generated during the PT rearrangement 
of the original rhs of Eq. ([2]) . 

A new SD equation for the gluon propagator - The relevant PT rearrangements take place 
when the longitudinal momenta of the three-gluon vertex trigger the STIs satisfied by specific 
subsets of fully dressed vertices appearing in the ordinary perturbative expansion. Unlike 
QED, due to the non-linearity of the BRST transformations, these STIs are realized through 
auxiliary (ghost) Green's functions involving composite operators such as (0|T[s$(x) ■ • • |0), 
where s is the BRST operator and $ is a generic QCD field. It turns out that the most 
efficient framework for dealing with these type of objects is the so-called Batalin-Vilkovisky 
formalism 12|. In this framework, one adds to the original gauge-invariant Lagrangian Ci 
the term £brst = J2-s> ^*s^, coupling the composite operators 5$ to the BRST invariant 
external sources (usually called anti- fields) $*, to obtain the new Lagrangian £bv = -^^i + 
-Cbrst- One advantage of this formulation is that it allows one to express the STIs of 
the theory in terms of auxiliary functions which can be constructed using a well-defined 
set of Feynman rules (derived from £brst)- In particular, the usual STI satisfied by the 
three-gluon vertex, an essential ingredient in the ensuing construction, assumes the form 



'?°TAjiA™A[j(^l, ^2) 



^£"'^1^717(^2, A;i)r^d^m(A;i)r^a'^™^*7(A;i, k2)TAdAn 



ih) , (3) 



where r^a^b(g) = (A ^)„^(g) - i^'^^qaqp [with -TA^Ap{q) = Ua/siq)], and the auxihary 
function Tcaa*, given in FiglJK, is nothing but the standard function appearing in the con- 
ventional derivation [l^ now written in the anti-field language. An important property of 
auxiliary functions involving the gluon anti-field. A*, is encoded into the so-called Faddeev- 
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--►--•=H = -S''"'k'y + 



4*7 



A7 



A?, 



^d 



-r'„ 



At 



At 



FIG. 2: The auxiliary functions — r^m^*7, — rQa^*7, c°-' a^a*~' ^"^^ ^^f^gA^g"- Black and white 
blobs represent one-particle irreducible and connected Green's functions, respectively, while gray 
blobs are connected kernels. 

Popov equation: ^ + iq^-^j^ = 0. This equation amounts to the simple statement that 
contracting A* with its own momentum q converts it to an anti-ghost, c. This property will 
be used extensively in what follows. In addition, one can obtain a set of useful identities re- 
lating Green's functions of background fields to those of quantum fields. These Background 
Quantum Identities (BQIs) llj are realized through auxiliary functions involving normal 
fields, anti-fields, and a background source fl, coupled through the term —gf'^"^"'c°'Q'^At^, 
see Figj2]D. The BQIs satisfied by the gluon propagator are 



iT 



4 a Ab 



which can be combined into the single identity 



^r^a^t — 'i^A'i.Al + ^^%A*;'^ AiAl + ^Q^AV^ A<XAi 
+ rf^a^*7i AiAl^ CM« 



/3 '"a^'-d ^l^fi ^'•f3^d '^a'^^ 

Y Ad AeV c 



Other BQIs needed in our construction will be 
where {ip, 0) G {{A, A), (c, c), (c. A*)}, and 



Rn^c^A 



— 4*7 Fr 



c^A*^-"- Q.%Ai^At^ 



(4) 
(5) 



(6) 



(7) 



(8) 



Equipped with these relations we may now proceed to the derivation of our main result. 
The aim will be to start from the conventional SDE of Fig{T^ and generate dynamically 
through the PT algorithm all the terms appearing in the BQI of Eq.([n]), thus arriving at the 
SDE equation of Fig{TJl. This will be accomplished by constructing the two BQIs of Eq.(jl]) 
and one at a time. The starting point is diagram (ai) of Fig{T^. The tree-level three 
gluon vertex F can be decomposed ^ into the sum F^ + F^, where (factoring out the color 
structure) 



'^^A^A^aA^I^ ^2) = gavkl^,- Qap^hv (9) 

This splitting assigns a special role to the physical momentum g, making F^ Bose symmetric 
only with respect to the and A^, legs inside the loop. In fact, F^ coincides with the BFM 
vertex F*^'' , , . F^ contains the longitudinal momenta that will get contracted with the full 

AAA ^ 

three-gluon vertex, triggering the STI of Eq.([3]). The result will be (ai) = (o-i)^ + (oi)^, 
with {aiY coinciding with diagram (61), and 

+ / D{ki)V'^,j^,^i{k2)D{k2)V^mAyA-(l^k2) +i I Z^(/i;i)F^™^6^,d(-g, /i;2) L (10) 

J ki J ki ) 

In the equation above we have used the ghost SDE k'^D{k) = 1 — iV cc{k)D{k) to transform 
a tree-level ghost propagator appearing in the second term of the rhs into a full one. The 
first integral on the rhs of Eq. lITO!) symmetrizes the ghost-gluon vertex of (03), giving rise 
to the characteristic BFM vertex oc {ki — ^2)0, and thus to diagram (63). The second term 
coincides precisely with the diagram (610); the third term (see Fig 12]) gives rise to diagram 
(64) [through the tree-level part of Fc^^*], as well as (fog), (^7) and (63). Finally, due to 
the fact that the four gluon vertices ^ Xaaa ^"^^ ^aaaa coincide at tree-level, we will have 
(0.2) = (^2), (0.4) = (^5), and (as) = (foe)- Thus taking into account the first term in Eq. (fTO|) 
we have dynamically reproduced the propagator BQI of Eq.(jl]). 

At this point we have constructed F^^^t (g); the next step will be to exploit the obvious 
equality F j„ .b iq) = ^ Ao-Ab (q) to interchange the background and quantum legs (see FigIT}:;). 
This introduces a considerable simplification: on the one hand we keep identifying the 
pinching momenta from the the PT decomposition of the (tree-level) F, while on the other 
hand the equality between diagrams (cs), (ce) and (cis), {do) is immediate. 



Let us now carry out the PT splitting of Eq.Q to diagram (ci). The part of the vertex 
generates directly diagram (rfi); the longitudinal momenta contained in get contracted 
as before with the full three-gluon vertex, which, however, has now an external background 
leg. Using Eq.([7]) with {(p,(f)) = (^4,^4), we get 

The presence of the prefactor ig'^d^'^ + T^J^a*^ allows one to use the BQIs of Eq.d?]) to convert 
the full vertices T^ac and Fcaa*, appearing in the last three terms of {a\), into F^^g and 
^cAA'^ respectively. This operation has two effects: {i) it generates (rfy), (dg), (dg), and 
(c/io), plus the contribution needed to convert (cs) and (C4) into (^3) and (^4), respectively; 
(m) it gives rise to leftover contributions given by the three integrals appearing in Eq. (fTO|) 
where the corresponding vertex is replaced by either —Rqcc or —RqcA*- These latter terms 
cancel exactly against the second term in Eq. (fTT!) . after its tree-level contribution has been 
extracted and used to convert (02) into (^2). At this point we have generated all diagrams 
of FiglDi In addition, using the BQI of Eq.(jl]) (already proven in the previous step), 
the term in {a\) proportional to Tqa* will give precisely —iTQaj^*i{q)T^dXb^{q)- Thus, we 
have constructed the full BQI of Eq.Q. Having dynamically realized the BQIs of Eqs.Q 
and ([5]), we can combine them into Eq.([6]), which constitutes the announced result. We 
emphasize that (z) all rearrangements have been induced by the PT manipulation of only one 
diagram [(ai) and (ci) of Fig{T] and (ii) all quantities encountered exist in the conventional 
formulation. In that sense, the Batalin-Vilkovisky formalism serves simply as an efficient 
way of keeping track of them. 

Discussion - The new SD series just constructed reads 
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[1 + GiqY^^M = (12) 

i=l 

where G is defined as the part of Fq^^* proportional to g^is- ^n^A*{q) = iG{q^)gap+ ~ qaQp- 
The rhs of Eq. (1121) has a very special structure. The diagrams of Fig|T]i can be separated into 
four subgroups [(di), (^2)], [(4), (4)], [(4), (4)], and [(dr), (cis), (c^g), (c^io)], corresponding 
to one- or two-loop dressed gluonic or ghost contributions. Due to the abelian Wis satisfied 
jy these new vertices, the contribution of each of the four subgroups is individually transverse 



15|. 
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The practical implications of this property for the treatment of the SD series are far- 
reaching, since it furnishes a systematic, manifestly gauge- invariant truncation scheme. In 
the case of the gluon self-energy, for instance, the transversality of the answer is guaranteed 
at every step. Specifically, keeping only the diagrams in the first group, we obtain the 
truncated SDE 

n„^(g) = [1 + G{q^)]-^[{di) + {d2)]a,3, (13) 

and we have that g"[((ii) + {d2)]ap = by virtue of q°T ^a^rn^nih, h) = - 
A~J(fc2)]- Therefore, Ilaf^{q) is transverse, as it should, despite the omission of the remaining 
graphs (most notably the ghost loops). In fact, one can envisage the possibility of employing 
completely different treatments for each subgroup: for example, one may treat the graphs 
[di) and (^2) non-perturbatively, while opting for a perturbative treatment of the ghost 
diagrams (^3) and (^4), without compromising the transversality of the self-energy. The 
price one has to pay is the need to consider the additional SDE governing G (see Figj2)3). 
Notice, however, that the approximations employed for the treatment of this latter SDE will 
not interfere with the transversality of Ilaf3- The abelian Wis furnish an additional technical 
advantage: one may use gauge-technique inspired Ansatze, a common practice when dealing 
with the SDE of QED to express the vertices in terms of propagators, in such a way as 
to automatically enforce gauge invariance. Finally, notice that (i) the SDEs for the QCD 



vertices can be constructed in a very similar way 16|], and (ii) the analysis presented here can 



be generalized to other gauges {e.g., the Landau gauge) using the methodology developed 
in Id]. 

In conclusion, the new SD series constructed in this letter provides a powerful tool for 
the systematic exploration of the non-perturbative sector of QCD, allowing the study of the 
fundamental Green's functions in a manifestly gauge- invariant way. 
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